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The Effect of Oscillating Fermi Energy on the Line Shape of the Shubnikov-de Haas 

Oscillation in a Two-Dimensional Electron Gas 

Akira Endo * and Yasuhiro Iye 

Institute for Solid State Physics, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8581 

The line shape of the Shubnikov-de Haas (SdH) oscillation has been analyzed in detail for 
a GaAs/AlGaAs two-dimensional electron gas. The line shape, or equivalently the behavior of 
the Fourier components, of the experimentally observed SdH oscillation is well reproduced by 
the sinusoidal density of states at the Fermi energy that oscillates with a magnetic field in a 
saw-tooth shape to keep the electron density constant. This suggests that the broadening of 
each Landau level by disorder is better described by a Gaussian than by a Lorentzian. 

KEYWORDS: Shubnikov-de Haas oscillation, two-dimensional electron gas, higher harmonics, Fermi en- 
ergy, Gaussian, Lorentzian 
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1. Introduction 

The Shubnikov-de Haas (SdH) oscillation is a preva- 
lently observed phenomenon in two-dimensional electron 
gases (2DEGs) subjected to a magnetic field perpendic- 
ular to the 2D plane. The oscillation is periodic in 1/B 
with the frequency proportional to the electron density 
Tie, thus serving as a standard tool to evaluate rif,.^ The 
damping of the oscillation with decreasing magnetic field 
reflects the scattering of electrons out of the cyclotron 
orbit, and therefore is widely used as a measure of the 
single particle scattering time tq or the corresponding 
quantum mobility /ig = erq/m*,^ with m* the electron 
effective mass. 

Despite the ubiquity of the SdH oscillation in varieties 
of magnetoresistance experiments on a 2DEG, practi- 
cally only the two aspects, the 1/B frequency and the 
damping, have been exploited so far to characterize the 
2DEG used in the experiment. This requires examining 
only the fundamental component of the oscillation, i.e., 
the sinusoidal oscillation with a due exponential damping 
factor. To the best of our knowledge, no report has been 
made to date on detailed investigation of the line shape, 
the behavior of higher harmonics,'^ of the SdH oscillation 
in a 2DEG.4 

The SdH oscillation derives from the oscillating den- 
sity of states (DOS) mainly through the modification of 
the scattering rate. The oscillation of the DOS, in turn, 
originates from the Landau quantization. 
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with LUc = eB/in* , which turns the energy independent 
DOS at _B = into a set of delta function peaks at 
E ^ En for an ideal 2DEG. The Landau level (LL) 
peaks, in reality, acquire width due to impurity scat- 
tering, which is usually modeled either by a Lorentzian 
P^{E) = (rL/7r)/(i;2 r^), or by a Gaussian Pg{E) = 
(l/%/2^rG)exp(-£;V2r^). The line shape of the resul- 
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will certainly be reflected in the line shape of the SdH 
oscillation. In eq. (2) we included the factor 2 to account 
for the spin degeneracy, a = L (G) for Lorentzian (Gaus- 
sian) peaks, and I = ^Jh/eB represents the magnetic 
length. The analysis of the line shape of the SdH oscil- 
lation, therefore, allows us, in principle, to gain insight 
into the LL peaks that constitute the DOS. 

Detailed knowledge of the DOS or the constituent LL 
peaks is indispensable to the quantitative understanding 
of a multitude of phenomena that originate from the Lan- 
dau quantization. The two models of the LL broadening, 
the Lorentzian and the Gaussian, considerably differ at 
their tails. Therefore their difference can be crucial in 
the quantitative interpretation of the phenomena that 
take place at LL tails (e.g., localization in the quantum 
Hall states), as well as of those that occur in the low- 
magnetic field region where adjacent LLs substantially 
overlap. A number of experimental techniques have been 
applied to the exploration of the DOS of a 2DEG, includ- 
ing the measurement of the specific heat,^ photolumines- 
cence,^ magnetocapacitance,^'® and magnetization. ^"^^ 
In comparison with these techniques, magnetoresistance 
measurement can be carried out with a simpler experi- 
mental setup thus generally with higher signal-to-noise 
(s/n) ratio, which is advantageous in investigating low 
magnetic-field range where the amplitude of the oscilla- 
tion is expected to be small. A major drawback of the 
magnetoresistance in this respect is the possible interven- 
tion by effects other than that of the DOS; these include 
the weak localization effect, localization in the quantum 
Hall regime, formation of the edge states. These difficul- 
ties can mostly be circumvented by limiting ourselves to 
low magnetic fields where adjacent LLs have sufficient 
overlap to prevent the localized states in the bulk, hence 
also the edge states, from being generated. Weak local- 
ization, if any, is already suppressed at the magnetic field 
at which the SdH oscillation appears in the modern high- 
mobihty GaAs/AlGaAs 2DEGs. 
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In the present paper, we make an analysis of the ex- 
perimentaUy observed hne shape of the low-field SdH os- 
cillations. Comparison is made with the line shape of 
calculated DOS. Emphasis is on the importance to take 
into consideration the oscillation of the Fermi energy Ep 
that keeps the electron density constant. 

2. Experimental Details 

We examined several GaAs/AlGaAs single- 
heterostructure 2DEG wafers with slightly varying 
sample parameters but all grown with the same molec- 
ular beam epitaxy (MBE) chamber. They all have 
a structure from the surface: 10 nm GaAs cap / 40 
nm Si-doped {^10^^ m~^) Alo.33Gao.67As / undoped 
Alo.33Gao.67As spacer having the thickness dg = 40 or 60 
nm / GaAs with 2DEG channel residing at the interface 
with the spacer layer. The electron density and the 
mobility ^i are 1.5-2.0x101^ m^^ and 50-70 mVVs, 
respectively, in the dark, which increase to 2.2—3.0 x 
10^^ m~^ and 70—130 m^/Vs after illumination by an 
infrared light emitting diode. The samples are fabricated 
into Hall bars for magnetotransport measurements. All 
the measured SdH traces taken from different samples or 
with different illumination conditions showed essentially 
the same features, with minor variations to be mentioned 
when necessary. In what follows, we present a typical 
example taken at the conditions ~ 2.9x10^^ m~^ 
and /i = 74 m^/Vs obtained after illumination. 

The measurement was carried out in a dilution fridge 
at the base temperature (~ 15 mK). We employed stan- 
dard low-frequency (13 Hz) ac lock-in technique for resis- 
tance measurement, with an excitation current of /rms = 
10 nA; the only effect found by reducing Ii-^s down to 0.5 
nA was to worsen the s/n ratio, attesting to the absence 
of the current heating. In order to capture the finest de- 
tails of the oscillation, we adopted a slow sweep rate (0.01 
T/min) of the magnetic field and a rapid data acquisi- 
tion rate (~ 4 data points/s), which amounts to roughly 
1 data point per 4xl0~^ T. The slow sweep rate is also 
beneficial in avoiding undesirable hysteresis of the su- 
perconducting magnet. The applied magnetic fields were 
further calibrated by the simultaneously measured Hall 
resistivity. 

3. Experimental Results 

Figure 1 shows a typical low-field magnetoresistance 
trace, with x the direction of the current and po the 
resistivity pxx aX B = 0. The amplitude of the SdH 
oscillation monotonically increase with increasing B up 
until ~0.3 T. The deviation from this trend at higher 
magnetic fields is attributable to the onset of the spin 
splitting. In the present paper, we examine only the low 
magnetic-field region where spin degeneracy remains un- 
resolved. Note that there the bare Zeeman splitting is 
only a few percent of the disorder broadening Fq of the 
LLs owing to the smallness of the g factor in GaAs,^^ 
and is therefore completely negligible. The spin splitting 
that evolves into the odd-integer quantum Hall states 
at higher magnetic fields is predominantly due to the ex- 
change enhancement of the g factor that takes place only 
above a certain magnetic field; the onset is considered to 
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Fig. 1. (Color online) A low-field magnetoresistance trace of a 
2DEG measured at ~ 15 mK. Upper/lower envelope curves and 
their average are shown by dotted lines and a dot-dashed line, 
respectively. The horizontal bar indicates the range of the mag- 
netic field plotted in Fig. 2. 
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Fig. 2. The SdH oscillation extracted from Fig. 1, plotted against 
the filling factor. The corresponding magnetic field is indicated 
by the top axis. 



be a second-order phase transition-'^^ and consequently 
the exchange enhancement is irrelevant below the tran- 
sition field. In Fig. 1, we also plot upper and lower enve- 
lope curves (dotted lines) and their average (a dot-dashed 
line). The envelope curves are obtained as spline curves 
connecting the extrema. The oscillatory part of the mag- 
netoresistance Apsdn/po is extracted from Pxx/po by 
subtracting the average curve. The resultant SdH oscil- 
lation ApsdH/[po^{T /Tc)] is plotted in Fig. 2 against 
the LL filling factor i' = n^h/eB for the magnetic-field 
range depicted by a horizontal bar in Fig. 1. The factor 
A{T/Tc) with fcsTc = (l/27r2)fiwc and A{x) = 2;/sinh(a;) 
corrects for thermal damping; the decrement of the fac- 
tor from unity is less than 0.01 for the relevant magnetic- 
field range, indicating that thermal damping is negligibly 
small for our low temperature. 

The SdH trace in Fig. 2 takes minima and maxima 
at even and odd filling factors, respectively, as is ex- 
pected for a spin-degenerate 2DEG. A notable feature 
to be highlighted in the present paper is the asymmetry 
between the maxima and the minima. While the peaks 
at the maxima exhibit rather dull rounded inverted U- 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 3 




10 15 20 

B ' (T ') 



Fig. 3. (Color online) Inset: the Fourier spectrum obtained from 
the Apg^ii/[poA{T/Tc)] vs. 1/B curve. Main panel: the ampli- 
tudes of the fundamental (red squares), the second (green upward 
triangles), and the third (blue downward triangles) harmonic 
contents of the SdH oscillation [see Fig. 4(a)] plotted against 
1/B. Amplitudes of the total SdH oscillation (Fig. 2) are also 
plotted by open squares. The line with the least steep slope rep- 
resents the fit of C exp(— 7r//tq_B) to the total SdH amplitude 
(after omitting the higher magnetic- field range), which defines 
^Q=lf.9 m^/Vs. The two steeper lines show C exp(— 27r//iQ_B) 
and C exp(— 37r//iQ_B) with the same values of C and fiQ as the 
first line. The horizontal bar indicates the magnetic-field range 
plotted in Figs. 2 and 4. 



shape, the dips at minima take on sharp V-shape. This 
trait in the hne shape is a quite generic feature observed 
in aU the 2DEGs we have investigated and also can be 
seen in the papers pubhshed by other authors (see, e.g.. 
Fig. 1 in Ref. 14). 

In order to quantify the hne shape, we carried out 
Fourier transform to the trace of ApsdH/[/Oo^(T/7c)] 
vs. 1/B (note that 1/B oc v so that the trace to be 
Fourier transformed is basicaUy the same as the trace 
shown in Fig. 2). The Fourier spectrum shown in the in- 
set of Fig. 3 exhibits peaks deriving from the fundamen- 
tal periodicity and up to the third (and a small trace of 
the fourth) harmonics. The resistivity component corre- 
sponding to each Fourier peak (up to the third harmonic) 
is obtained by performing a Fourier band pass filter to 
the ApsdH/[po^(J"/Tc)] vs. 1/B curve using the shaded 
region in the Fourier spectrum as a window, and plotted 
in Fig. 4(a) against ly. Addition of the three traces in Fig. 
4(a) practically reproduces the trace in Fig. 2. 

The amplitudes (absolute values of maxima/minima) 
of each oscillatory components in Fig. 4(a) are plotted 
in the main panel of Fig. 3 (solid symbols) in the semi- 
logarithmic scale, along with the amplitudes obtained 
from Fig. 2 (open squares). It is clear from the figure 
that the SdH oscillation is dominated by the fundamen- 
tal component. This justifies the conventional treatment 
in which the oscillation is approximated by a single sinu- 
soidal curve with an exponential damping factor 
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The quantum mobility /iq is deduced from the damp- 
ing of the amplitude. We obtain /iq = 11.9 m^/Vs 
by the fit of C exp{—7T / iiqB) to the total SdH ampli- 
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Fig. 4. (Color online) (a) The fundamental (red), the sec- 
ond (green), and the third (blue) harmonic content of the 
SdH oscillation (plotted with progressively thicker line) ob- 
tained by a numerical Fourier band pass filter performed on 
Aps^ll/[poA{T/Tc)] shown in Fig. 2. The frequency windows 
allowed by the filter are indicated by the shade in the inset of 
Fig. 3. (b)(c) Similar to (a) with the band pass filter performed 
on the Lorentzian DOS at fixed Ep [Fig. 5(a)] (b), or on the 
sinusoidal DOS at oscillating Ep [Fig. 6(b)] (c). Note the dif- 
ference in the phase of the second harmonic component between 
(b) and (a), (c). See text for detail. 



tudes, which fits to the amplitudes of the fundamen- 
tal component as well, as seen in Fig. 3. The higher 
magnetic-field regime where spin-splitting commences is 
omitted from the fitting. Interestingly, the amplitudes 
of the second and the third harmonics roughly fall on 
the fines Cexp(— 27r/^qi3) and Cexp(— 37r//iqi?), re- 
spectively, with the identical values of C and fiQ ob- 
tained by the above fitting. This property is shared by 
the SdH traces taken at other illumination conditions 
(i.e., different and rig) and with other samples. "'^^ As 
will be discussed later, this behavior of the amplitude 
is what we expect for a Lorentzian line shape, although 
the Lorentzian fails to reproduce the correct phase of the 
second harmonic. In the subsequent section, we compare 
our experimental line shape with two types of calculated 
DOS. 
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4. Comparison with Calculated DOS 

4..1 Oscillation of the Fermi energy 

For low magnetic fields, it is convenient to rewrite, 
with the aid of the Poisson sum formula, the DOS 
eq. (2) in the Fourier series, 

D{E) = 



Do\l + 2j2 



cos 



fe=l 



E 



X exp 



(4) 



with Do = m* jixT? the constant DOS at i? = 0, and 
p = 1 for Lorentzian and 2 for Gaussian LL peaks. Here 
we resorted to an approximation Y^=q ~^ X^?/=-oo' 
approximation causes virtually no change in the DOS at 

^ Fq,, i.e., at E beyond the reach of the tail from the 
peaks located at A'' < 0. Since terms with larger fc decay 
more rapidly owing to the exponential damping factor, it 
is usually enough to take only a few terms into account 
in the summation of eq. (4) for low magnetic fields. 

Most of the LL peaks experimentally measured so far 
have been explained either by a Lorentzian with the E- 
independent width''' or by a Gaussian with the width 
proportional to \fB^ With 

rL = ;^- (5) 



2to* 



and 



eh 
2m* 



/ 2B 
tt/xq' 



(6) 



to 



the exponential factor in eq. (4) reduces 
exp(— fc^'Tr/yUqi?), and by substituting k = 1 coin- 
cides with the exponential factor in eq. (3) for both 
Lorentzian and Gaussian. 

Using the D{E) in eq. (4), the electron density at low 
temperatures is written as 

rie = / D{E)f{E)dE ~ / D{E)dE = N{Ey) (7) 
Jo Jo 

with f{E) = {l + cxp[(E-EF)/kBT]}~^ the Fcrmi-Dirac 
distribution function. We defined the cumulative number 
of states (i.e., integrated DOS) below an energy E as 



N{E) 
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.. (8) 



The electron density of a 2DEG is expected not to 
vary with magnetic field at low temperatures. To keep rig 
constant, eq. (7) requires Ep to oscillate with B around 
the value at B = 0, Ep = Ue/DQ. The oscillation of 
Ep with B, as well as the DOS at the oscillating Ep, 
is evaluated for the two types of LL broadening in the 
following subsections. 
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Fig. 5. The upper two panels show the oscillatory part of the 
Lorentzian DOS, with the width Fl determined from the exper- 
imentally obtained /iq =11.9 m'^/Vs. (a) assumes a fixed E^. = 
10.4 meV (corresponding to the experimental Ue = 2.92x10^' 
m^^), while (b) takes account of the oscillating Ep calculated 
by the procedure described in the text, (c) shows asymmetric 
(sawtooth-like) oscillation of the Fermi energy, repeating alter- 
natingly steeper rises and gradual declines with i/. 



4-2 Lorentzian density of states 

The oscillatory part AD /Do ^ (D - Do) /Do of the 
Lorentzian DOS at fixed Fermi energy E^ is plotted in 
Fig. 5(a). The line shape with sharp peaks and rounded 
dips is at obvious variance with the experimental line 
shape of the SdH oscillation. Equation (4) along with 
eq. (5) shows that the amplitude of the A;-th harmonic 
is given by 2 exp(— fcTr/^qi?). Therefore the decay with 
1/B of the amplitude of the harmonics for the experi- 
mental SdH oscillation presented in the previous section 
(see Fig. 3) is in accord with that of the Lorentzian. How- 
ever, the oscillation of each harmonic content shown in 
Fig. 4(b) obtained by the Fourier band pass filter reveals 
that the phase of the second harmonic of the Lorentzian 
DOS is inverted from that of the experimental trace. 
Therefore, the experimental SdH oscillation is more like 
a Lorentzian laid upside-down, AD —AD after a shift 
Ep ^ Ep + hLOc/2 [see eq. (4)]; it looks as if the gaps be- 
tween LL peaks are comprised of inverted Lorentzians. 

As discussed in the previous subsection, it is necessary 
to take the oscillation of the Fermi energy into account. 
For the Lorentzian DOS, AEp can be calculated analyt- 
ically. Equation (4) with p = 1 can further be rewritten 
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as 



Die) = Do- 



sinh(27r7) 



(9) 



cosh(27r7) — cos(27re) 

Here we introduced dimensionless parameters, 7 = 
Tl/ZiWc = 7r(rG/?iWc)^ = (2mq-B)-^ and e = E/HuJc - 
1/2. In what follows, we will also use dimensionless 
Fermi energies ep and £p similarly defined as e, and 
also d = AEp/fuJc for the deviation AEp = Ep - 
of the Fermi energy from the zero-field value. Note that 
£p = 1//2 — 1/2. Accordingly, the cumulative number of 
states reads 

N{e) = 

2^ i — arctan [coth(7r7) tan(7re)] + i + Int 

h [ TT 2 



r 

^+2 



where Int[a;] signifies the integer part of x. From eqs. (7) 

and (10), one obtains 

tan(7r£p) = coth(7r7) tan(7r£F) 

= coth(7r7) tan [71(4 + S)] , (11) 
leading to the oscillatory part of the Fermi energy 
1 



S = - arctan | [tanh(7r7) - 1] tan(7r4) | _ 
TT I 1 + tanh(7r7) tan^(7rep) 



(12) 




The resulting AEp is plotted in Fig. 5(c). Prom eq. (11) 
one obtains 

. . 1 — tan^(7r£F) 

C0S(27r£F) ^ 



1 +tan2(7r£F) 
1 - tanh^(7r7) tan2(7r£^) 



1 + tanh^(7r7) tan^iirs^] 

The DOS at oscillating Ep is obtained by replacing eq. 
(13) into eq. (9) 



D(£f) = Do 



coth(27r7) + 



1 



■ cos(27r£^) 



(14) 



sinh(27r7) 

whose oscillatory part is plotted in Fig. 5(b). Interest- 
ingly, higher harmonics arc gone and the oscillation in- 
cludes only the fundamental component. To conclude 
this subsection, the Lorentzian DOS does not reproduce 
the line shape of the experimentally observed SdH oscil- 
lation, even with the oscillatory Ep. 

4.-3 Sinusoidal density oj states 

For the Gaussian DOS, it is rather difficult to perform 
similar calculations as were done in the previous subsec- 
tion for the Lorentzian DOS. Instead, we study the sinu- 
soidal DOS in this subsection. Since the terms with k > 
2 in eq. (4) decay rapidly for the Gaussian DOS because 
of the factor fc^ in the index of the exponential damping 
factor, the sinusoidal DOS that keeps only the term k = 
1 in the summation of eq. (4) constitutes a good approxi- 
mation for the Gaussian DOS. Actually, the amplitude of 
the term k = 2 accounts for only 0.1— 3.8 % of that of the 
term k = 1 for the current example in the magnetic-field 
range of the present interest. For a sinusoidal DOS 

D{e) = Do [1 -I- 2 cos(27r£) exp(-27r7)] (15) 



Fig. 6. Similar to Fig. 5 for the sinusoidal DOS, which serves as 
a good approximation for the Gaussian DOS as explained in the 
text. 



(13) 

we obtain 



N{e) = 2 



eB 



£ -h ^ -h — sin(27r£) exp(— 2777) 

2 TT 



Equations (7) and (16) result in 

= ,5 + - sin [2tt {e'^ + 6)] exp{-2TTj). 



(16) 



(17) 



Equation (17) may be solved numerically. Here, instead, 
we deduce an approximate solution valid up to 0{6^), 
noting the smallness of S in the magnetic-field range of 
the current interest. We get 



5 ~ 



1 + 2Xc° - V(l + 2AcO)2 + 2(2AsO)2 



(18) 



47rAsO 

where we have introduced notations A = exp(— 2777), 
c° = cos(27r£p), and s° = sin(27r£p) for brevity. The 
calculated AEp is plotted in Fig. 6(c) along with the os- 
cillatory part of the DOS either at fixed Ep or at oscillat- 
ing Ef in Figs. 6(a) and 6(b), respectively.^^ Figure 6(b) 
reproduces the line shape of the experimental SdH oscil- 
lation in Fig. 2 quite well. Both the amplitude and the 
phase of the harmonics are basically reproduced, as can 
be discerned by comparing Figs. 4(a) and 4(c), the latter 
showing the harmonic contents of Fig. 6(b) obtained by 
the Fourier band pass filter. 

To be more quantitative, we deduce approximate for- 
mula for AD /Do at oscillating Ep upto 0{X'^). Note that 
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the sinusoidal DOS starts to deviate from the Gaussian 
DOS only at 0(A''). Substituting 

2tt5 ~ -2sin(27r4)A + 2 sin(47r4)A^ + 0{\^) 



into 



^^^^^^ = 2 cos(27r£F)A = 2 cos [27r(4 + 5)1 A, (19) 



one obtains 

A£)(£F 

Do 



= 2 cos(27reSi) A + 2 [l - cos(47reS^)] A^ 

+3 [- cos(27r4) + cos(67r4)] A^ + 0{X^). 

(20) 

Equation (20) gives the right amplitude and phase of the 
harmonics as inferred from Figs. 3 and 4(a), except for 
the 1.5 times larger amplitude of the third harmonic. The 
behavior of high harmonics, however, is rather subtle. In 
fact, the use of eq. (12) instead of eq. (18), a fairly good 
approximation considering the minuteness of the differ- 
ence between Figs. 5(c) and 6(c), leads to the expected 
amplitude and phase up to the third harmonic (see Ap- 
pendix). 

5. Discussion 

Different experimental approaches employed to ex- 
plore the LL line shape^"^^ differ in the preference 
of Lorentzian or Gaussian. Recent experiments,'^' 
however, seem to converge on the Lorentzian broadening 
with Fl independent of B for middle to high magnetic 
fields, B >~1 T. For lower magnetic fields, the result of 
recent high sensitivity magnetization measurement car- 
ried out down to a magnetic field as low as ~0.5 T can be 
fitted by the two types of broadening equally well, hence 
cannot distinguish between the two.^'^ Therefore there 
has been no general agreement on the LL line shape at 
low magnetic fields. The present study suggests that the 
Gaussian is the better candidate for still lower magnetic 
fields. 

Difficulty in magnetization experiments at low mag- 
netic fields appears to be arising at least partly from the 
limit in the sensitivity that hampers the acquisition of 
the data with sufficient s/n ratio. In magnetoresistance 
measurements, small amplitude oscillations can be de- 
tected with satisfactory s/n ratio as we have shown in the 
present study. On the other hand, the interpretation of 
magnetoresistance data becomes complicated at higher 
magnetic fields, affected by localization or edge states in 
the quantum Hall regime. Therefore the two experimen- 
tal techniques are more or less complementary. 

Theoretically, a pioneering work by Ando and Ue- 
mura^^ suggested semi-elliptic LLs with the width pro- 
portional to by self-consistent Born approximation 
(SCBA). The approximation is insufficient at the tail 
of LLs where multiple scattering plays important role, 
and the semi-elliptic LLs are not applicable at low mag- 
netic fields where adjacent LLs overlap. Alternative ap- 
proaches showed that the LL broadening is described by 
a Gaussian.^"' In these early studies, a short-range ran- 
dom potential is assumed. A later work^'' showed that a 



Gaussian line shape with the width Fq oc \/B holds also 
for a long-range random potential more appropriate for 
GaAs/AlGaAs 2DEGs. Our result is, therefore, in accor- 
dance with the theoretical prediction. 

In the comparison between the experimental SdH os- 
cillation and the calculated DOS, we have tacitly as- 
sumed the proportionahty A/9sdH/[po^] ^ AD/Dq- It 
has been claimed, with ample experimental evidence, 
that is proportional to the square of the DOS.^^'^^ 
This leads to the same proportionality as long as the 
oscillatory part AD of the DOS is much smaller than 
the constant background Dq so that an approximation 
{D/D^f ~ l + 2AD/Do is allowed. However, AD/Dq is 
not necessarily small enough even in the low magnetic- 
field range examined in the present study, as can be seen 
in Figs. 5 and 6. Therefore the relation p^x suggests 
that the comparison should be made between the exper- 
imental psdii/[poA] and the oscillatory part of D'^. We 
have actually made such comparison, only to find much 
worse agreement with the experimental SdH oscillation 
for any combinations of the type of DOS and Ep (either 
fixed or oscillating). The relation p^x oc derives from 
the theory for short-range random potential. The ex- 
tension of the theory to the system with long-range ran- 
dom potential is not straightforward, involving subtleties 
in the treatment of relevant scattering times. Our re- 
sult rather suggests the relation ApsdH/[po^] oc AD/Dq 
remains valid regardless of the magnitude of AD/Dq. 
However, the possibility that this trait is specific to the 
samples we examined cannot be completely ruled out, 
since they are all grown in the same MBE chamber and 
therefore can possibly contain an unidentified common 
source of scattering unintentionally introduced during 
the growth. 

6. Conclusions 

We have made detailed Fourier analysis of the line 
shape of the SdH oscillation at low magnetic-field regime 
where the oscillation predominantly reflects the oscilla- 
tion in the spin-degenerate DOS. The line shape can for- 
mally be described by inverted-Lorentzian gaps between 
adjacent LLs, and can be reproduced by the sinusoidal 
DOS at the oscillating in a saw-tooth shape. The 
sinusoidal DOS is consistent with the Gaussian broad- 
ening of LLs with the width oc y/B, in agreement with 
the broadening theoretically calculated for high LLs in a 
smooth random potential. 

The analytic formulae for oscillating Ep that keeps the 
He constant are presented in eqs. (12) and (18) for the 
Lorentzian DOS and the sinusoidal DOS, respectively. 
The use of the oscillating Ep has been found to be crucial 
for the agreement between the line shapes of calculated 
DOS and the experimental SdH oscillation. 
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Appendix: Evaluation of Eq. (19) with Eq. (12) 

The saw-tooth variation of Ep rcsuhs from the grad- 
ual decrease in Ep with decreasing B (i.e., increasing 
v) while a (spin-degenerate pair of) LL is being filled, 
alternating with the sudden jump up to the next LL at 
(even) filling factors where the LL that has been the host 
of the Ep is filled up. The line shape of AEp, particu- 
larly the amplitude (the jump), is primarily determined 
by the Landau fan diagram eq. (1) and the width of each 
LL, and the detail of the LL line shape plays only mi- 
nor role of slightly altering the shape of the saw teeth. 
This explains the resemblance between Figs. 5(c) and 
6(c). Hence it seems to be an acceptable approximation 
to use AEp derived from another type of the DOS. Equa- 
tion (12) has an advantage of being written analytically 
without approximation. In the following, we will show 
that eq. (19) can readily be expanded up to an arbitrary 
order in A = exp(— 2777) by employing eq. (12) instead 
of eq. (18). The result lends itself to provide a simple 
unified view on the role of the oscillating AEp. 

First, we rewrite eq. (13) by using the relation 
tanh(7r7) = (1-A)/(1 + A), 

cos(27reF) = 
1 



A 



1 



l + Aco.s(27r4) 
l-|-A2-|-2Acos(27r4) 



+ A 



1 + Acos(27Tc;!0 



l + A2 + 2Acos(27r£' 



By employing an identity 



En , n\ 1 - a cos ( 
a" cos(n^) = 7; 
l + a2 

ra=0 



+ a? — 2a cos 9 ' 



we obtain 

cos(27r£F) = cos(27r£p) 

00 

+ E (-1)'' {cos [27r(fc + 1)e%] - cos [27r(fc - 1)4] } A'= 

fe=i 

(A 

to be replaced in eq. (19). If we pick out the main con- 
tribution (the lowest order in A) for each harmonic, we 
get 



ADjep) 

Dn 



-2 ^(-1)'= cos (27rA;4) A'^ (A-3) 



fe=i 



or 



DiEp) . Do |l - 2f:cos (2./cg) exp (-2^I^fc) | 

(A;4) 

which is nothing but a shifted and inverted Lorentzian 
DOS (see the discussion in the first paragraph of §4.2). 
Therefore, to summarize, the oscillating AEp turns the 



Lorentzian DOS into a sinusoidal line shape, as described 
in §4.2 [eq. (14)], and, in turn, transforms the sinusoidal 
DOS into an inverted Lorentzian line shape. 
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